Numerical Studies of Quantum Oscillations in the Superconducting vortex mixed state 
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(Dated: December 17, 2008) 

We studied the quantum oscillations in superconducting vortex mixed states with d-wave pairing. 
We showed that the Onsager relation does not always apply. Furthermore, at mean field level the 
quantum oscillations from electron pockets are suppressed by the pairing. We conclude that an 
interpretation ol the experimental results asccoming from the four hole pockets created by a (n, n) 
folding cannot be ruled out. 
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Recent experiments have shown quantum oscilla- 
tions in underdoped YBCO samples in strong magnetic 
Md~45T flSSSH. This observation has been inter- 
preted as a signal of the underlying Fermi surface (FS). 
The period of the oscillation is large, which corresponds 
to a very small FS, and implies that the FS has been 
reconstructed, probably by some translational symmetry 
breaking order in the ground state 0, 0, However, 
for both materials that have been studied, the simplest 
construction, a (-7T, 7r) folding forming four hole pockets 
in the original Brillouin Zone (BZ), would imply a pocket 
too small comparing to the nominal doping by about 25% 
P, [H . This partly motivated a number of workers to in- 
tepret the data in terms of more complicated reconstruc- 
tion, such as incommensurate spin density wave (SDW) 
@> 3- Furthermore, the measured Hall effect is nega- 
tive [10( and this has led LeBoeuf et al to propose that 
the quantum oscillations originate from electron pockets. 
Whether the negative Hall effect is due to flux flow 11 1 
is currently a subject under debate. 

We believe that in these experiments the samples are 
still in a vortex mixed state. One evidence is the mea- 
sured torque hysteresis [H, which implies that vortices 
exist at least up to 45 T. Furthermore, the commonly 
quoted core size of 20 A° (see for example, extrapola- 
tion based on STM measurements [13] as well as Nernst 
measurements [13]) lead to an estimate of of 100 T. 
Therefore, in order to interpret the oscillations data one 
has to understand the quantum oscillations in the mixed 
state. 

Up to now all discussions assume that quantum oscil- 
lations in the mixed state maintain the same frequency 
as in the normal state. This is true of all experiments 
performed up to date where it is possible to scan the 
magnetic field across H C 2 [3, 15 1. However, we note 



that all previous experiments have been done on conven- 
tional s-wave low T c superconductors, with the possible 
exception of the organics which may be d-wave, and the 



high T c cuprates may be in quite a different parameter 
regime. For example, the coherence length £o is very 
short, of order 4 or 5 lattice constants and is the conse- 
quence of a large energy gap Ao. The number of Landau 
levels N is about 10 in the high T c experiments, as op- 
posed to 100's or 1000's for coventional superconductors. 
This means that the semiclassical orbit encloses 10 flux 
quanta, i.e., 20 vortex cores. There are good reasons to 
believe that the pairing amplitude and the gap scale Ao 
are very robust in the underdoped cuprates. In contrast, 
the commonly used theory of Maki [l(| and Stephen 
considers the lowest order scattering of the normal state 
quasiparticles by a pairing order parameter which is ran- 
dom with short range correlations. It is not at all clear 
that this picture applies to the current problem. Another 
parameter which clearly distinguishes the cuprate from 
other superconductors is fc^Co which is less than or of 
order unity if Uf is the size of the observed pocket while 
^f£o >> 1 for conventional superconductors. 

In this paper we address the question of whether 
the traditional picture continues to hold in a parame- 
ter regime which has not been tested experimentally. We 
take as our model the Bogoliubov-de Gennes (BdG) equa- 
tions with a variety of vortex coordinates and competing 
order parameters, which can take on arbitrary spatial 
dependence. We take the coherence length and the pair 
field as parameters and make no attempt to solve the 
problem self-consistently. Since we do not expect the 
BdG equations to be the correct microscopic theory for 
the underdoped cuprates, it makes little sense to deter- 
mine these parameters self-consistently. Rather, we treat 
this as a phenomenological model. To the extent that the 
proximity to the Mott transition is not captured by by 
this phenomenology, the application of our theory to high 
T c problems should be treated with caution. With these 
caveates, we write down a tight-binding Hamiltonian on 
a square lattice constant a, and we set e — H = a = 1: 
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where t and t' are the hoppings, Ad is the nearest neigh- 
bor (NN) pairing, A s f is related to the staggered flux or 
equivalents d-density wave order (SF-DDW) [H, 
and V s , V c are the staggered spin and charge poten- 
tial respectively. The potentials are defined on sites 
fi = (ix,iy) and the pairing is defined on bonds f\j — 
(fi +r"j)/2. The ^-dependent d-wave gap would be of size 
2 Ad (cos k x — cos k y ), and near half filling the gap opening 
seen in the tunneling expriments would be ~ 4Ad- As a 
phenomenological model, we take t to be of order J, the 
antiferromagnetic coupling. In underdoped cuprates the 
gap size is then ~ 0.5t. Aij is the electromagnetic gauge 
field taken to be compact on the lattice, which satisfies 

Tplaqutte A H = B and Etriangle A V = B / 2 : B is the 

magnetic field. 

The pairing amplitude near a vortex is described by 
the ansatz 

|A d (r)| = A o sin0; 

y/r$ + d(r) 2 

where Ao is the NN pairing amplitude deep inside the 
superconductor, ro the core size, and d(r) is the distance 
to the vortex center. When multiple vortices are near by 
we replace d by d m i n = (J^i d\~ p ^)~( 1 / p ') which is smaller 
than the distance toward the nearest vortex. The choice 
of p > 1 does not qualitatively affect the result. The 
phase (j) of Ad is not a gauge invariant quantity so we 
must determine it according to the gauge choice of Ay. 
An alternative way is to assign (f> and then determine the 
gauge field under the constraint mentioned above. We 
choose <j> to follow the constraints: (i) |A0| < ir/2 for 
every link; (ii) ^2 loop A0 = ^ 7rn wnere n is the num- 
ber of vortices enclosed, and (iii) <fi is periodic in the y 
direction. We then determine Ay using the fact that 
the physical configuration must minimize the free energy 

^Aj(f)i) s (!f) 2 where v s — fa — fa — is the super- 
fluid velocity. The free energy is a function quadratic in 
A^, so we can optimize it by solving a linear equation, 
using sparse matrix routines. 

V c , A s f , and V s describe the order in the normal state. 
They can be uniform, periodic, or localized around the 
vortices depending on the order present. There is an ad- 
ditional piece in V c which balances the charge density in 
the mixed state. In a finite system with periodic bound- 
ary conditions it is not always possible to fit in a periodic 
vortex lattice. Instead we work with with a disordered 
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array of vortices. At each B the vortex positions are 
determined by a Monte Carlo annealing process for par- 
ticles with Cr -2 repulsion. The vortices are stuck as we 
gradually lower the temperature. The resulting vortex 
configuration has short range order but is rather disor- 
dered. It is a reasonable representation of a snapshot of 
a vortex liquid or a pinned vortex solid. 

In order to see the quantum oscillations, either the 
sample size has to be larger than the cyclotron radius, or 
one can take advantage of the periodic boundary condi- 
tions. With the periodic boundary conditions the mag- 
netic flux through the sample is quantized, and to achieve 
a small stepping in 1/B the sample size again has to so 
large that the hamiltonian cannot be diagonalized. We 
can instead, use an iterative green's function's method to 
get the local density states (LDOS), at any fixed energy: 

We affix our sample to two semi-infinite stripes in the 
±x direction. The stripes are normal metal described by 
t and t' . We take periodic boundary conditions in the y 
direction. Now the configuration is similar to [20( and we 
can use the same method: 

G L {x)= [G^x)- 1 -tG L (x- l)tt]- x ; (3) 
G(x) = 

(G^x)- 1 -tG L {x-l)tf -tfG R (x + l)t)-\ (4) 

where G (x) is the Green's function for the isolated Xth 
column, G L (x) (G R (x)) is the Green's function at col- 
umn x when the right (left) side of the column is deleted, 
t is the hopping matrix between the two consecutive 
columns (contains t,t' , A^{x), and A s f(ir)), and G(x) is 
the Green's function for the Xth column in the original 
setting. Everything here is a 2y x 2y matrix, containing 
the electron part and the hole part. We first go from left 
to right, compute G L for every x. Next we go from right 
to left to compute G R and use Eq ((U to compute G(x) 
for all x. The imaginary part of the yth diagonal matrix 
element in the electron part of G(x) is then related to 
the local density of states at (x, y). 

Using this method we can get the LDOS everywhere 
on our sample. However, the LDOS varies from place to 
place due to its relative position to the vortices as well 
as the disorder of the vortex lattice. To see the quantum 
oscillations, it suffices to look at the averaged density of 
states. In this work, if not specified otherwise we set 
t' = -0.3t, r = 5 and the lattice is of size 2000 x 80. We 
vary the magnetic field in accord with the number of the 



vortices, which is from 500 to 10000. 




FIG. 1: this plot shows the averaged DOS vs. 1/B in the 
vortex mixed state with (tt, tt) SF-DDW order. B is in units of 
( 2 ^° 2 ). If we view H C 2 as the field strength where vortices are 
overlapping, it corresponds to 1/B = 25 in the plot. Different 
lines in the plot show states with different pairing fields. The 
maximum pairing gap in the antinodal direction would have 
size ~ 4 A. We can see as the gap increases the frequency 
decreases. The inset is a plot of the Fermi surface in the 
normal state. There are 4 hole pockets with 2.5% area each 
of the original BZ. This corresponds to p — 0.1 from half 
filling. In this plot A sf = 0.25*. 

In Fig. [T] we showed the result where we start from 
the state in which the Fermi surface is reconstructed by 
(tt, 7r) SF-DDW order. If we define H c2 to be the field 
where vortices start overlapping and ro = 5a, this gives 
(sf^sOlfl = 25, where <j> here is the full flux quantum, 
hc/e. The experimental probe, roughly at the 10th Lau- 
dau level, would be around (^\)-ft— = 60. We create 

' V 2tra z * H c 2 

a gap large enough to kill the electron pockets, leaving 
4 hole pockets shown in the inset. We found that the 
oscillation period of the normal state matches the pre- 
diction from the Onsager relation and the Luttinger the- 
orem. As we turn on the d-wave pairing amplitude, the 
period of the quantum oscillation increases. As the pair- 
ing gets large the period of the oscillation is off from the 
period in the normal state by about 20%. The frequency 
shifts are clearly seen in the Fourier transform shown in 
Fig. [S] Note that in our modelling of the vortex cores, 
the superconducting pairing remains, albeit at reduced 
amplitude, even for B > H C 2. This explains why the fre- 
quency shift persists somewhat above H C 2, but we can 
ignore that region becouse it is not reached experimen- 
tally. We have also done calculations for a (tt, tt) spin 
density wave (SDW) order, and the period is shifted in a 
similar way (see Fig. |6(b)| ) 

To check whether this is a generic phenomena, we ran a 
simpler setting with an electron pocket centered at origin. 
The result is showed in Fig. [2] Again, as we turned on the 
d-wave pairing, the frequency is reduced as the pairing is 
increased. It is worth noting that we shifted the chemical 




FIG. 2: an electron pocket is centered at the origin. The 
pocket has an area of 14% of the original BZ. In this plot 
t' — — 0.14t. The maximum gap ~ 1.5A. 



potential in order to maintain the total electron density 
as we increase the pairing. We also added a charge po- 
tential to make the charge density approximately uniform 
inside and outside the vortex core. Even if we just forget 
about density and keep the same chemical potential as in 
the normal state, the period still decreases (not shown). 
The chemical potential shift is significant in the case of 
Fig. [2] because the pairing is strong and the electron is 
far from particle-hole symmtries. It is less significant in 
the more realistic case of Fig. Q] 
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FIG. 3: 4 hole pockets and 2 electron pockets are present. 
In this plot A h = 2.9% and A e = 1.44% of the original BZ, 
and p = 0.086. In the normal state, from Onsager's relation 
we can determine that the peaks with period 11 are from 
electron pockets. The amplitude is heavily suppressed in the 
mixed state. In this plot V s = 0.2t. 

In Fig. [3] we started with a 2-pocket Fermi surface 
reconstructed by a (tt, tt) SDW order, and turned on a 
small d-wave pairing. The oscillations from the electron 
pockets are clearly visible in the normal state, but are 
rapidly killed in the mixed state. Note the pair field is 
very small, about 10 times smaller than that in Fig. [T] 
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One explanation is that the electron pocket originates 
from the antinodal region where the gap is large, and 
is dephased by the random pairing potential due to the 
random vortex configurations. However, this cannot be 
the whole story because we find that for s-wave pairing 
with similar pairing gap size the electron and hole pockets 
both survive. One difference we noticed is that in the s- 
wave case there is a strong density of states peak inside 
the vortex core, much larger than that inside the d-wave 
core. At this point we do not have a full understanding 
of the suppression of the electron pocket oscillations. 

We have also calculated the oscillation for an "incom- 
mensurate" SDW, with period Q = ± 26), it) where 
6 = 1/8. We impose a sinusoidal V s (r) = V s q cos(2ir6x)) 
and V c (r) = V c o cos(4&r). Complicated band structures 
in this kind of potential was computed by Millis and Nor- 
man [7|. It is instructive to consider the hybridization 
only with the primary vector Q for SDW and 2Q for the 
charge component. Fig. [4] shows the pure SDW case 




FIG. 4: Primary bands for the incommensurate SDW in the 
first quadrant of the original BZ. 

(V c = 0) and we see that after hybridization only two 
closed orbits are possible. One is a small hole pocket and 
the other is the electron pocket. The larger hole pocket 
becomes an open orbit. Inclusion of V c further cuts the 
small hole pocket to an even smaller area. Fig. O shows 
the quantum oscillations in a uniform 1/8 stripe phase. 
We indeed see the oscillations from the electron pock- 
ets and the small hole pockets. The electron pocket is 
again heavily suppressed in the mixed state, leaving a 
hole frequency which is much too small compared with 
experiment. 

In conclusion we have shown that the Onsager rela- 
tion does not always apply to quantum oscillations in 
the vortex mixed state. When the pairing is strong and 
the coherence length short, we find a systematic decrease 
in the frequency. We have also checked that if the core 
size is increased, the shift is diminished. Thus our result 
is consistent with experiments performed so far on con- 
ventional superconductors, which are in the large core 
size limit. The implication for the experiment on high 
T c cuprate is that the simple interpretation in terms of 
(tt, 7r) folding creating four hole pockets cannot be auto- 
matically ruled out. There are two candidates for {%, n) 
order. One is SDW order and the second is SF-DDW 
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FIG. 5: 1/8 stripe state with t' = -OAt, V s0 = 0.4, and 
Vco = 0.1. In the normal state 3 frequencies are visible. The 
fast oscillations localized around 1/B ~ 20 is a breakdown 
effect. At lower fields the fast oscillation (period 6) is coming 
from the electron pockets and the very slow oscillation (only 
1 period seen in the range) is from the hole pockets. As we 
turn on superconductivity the oscillation from the electrons 
is diminished. 
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FIG. 6: Fourier transforms of plots shown in Fig[T][2] and [3] 

We take a = 4 A° and the frequency is then in units of Tesla. 
The color used is the same as in the averaged DOS vs. 1/B 
plots. The normal state frequency is shown by the green line, 
(a) SF-DDW.(b) SDW, p=0.125. The peak at the right is the 
second harmonic, (c) electron pocket centered at origin, (d) 
two-pocket SDW; the peak at ~380 T is from the electron 
pockets and the one at ~750T is from the hole pockets. 



order. We note that in view of the recent measurement 



that the g factor is much less than 2 [2l|], the SF-DDW 
scenario must be accompanied by some additional order, 
such as incommensurate SDW. In that case there are pos- 
sibilities for larger orbits composed of two or three hole 
pockets shifted by 8 if magnetic breakdown is taken into 
account. The other interpretation is the incommensu- 
rate SDW order. In this case the observed pocket must 
be identified with the electron pocket. Alternatives in 
terms of SF-DDW with electron porckets have also been 
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proposed. @, H, 0| . While the electron pocket is rapidly 
suppressed in the mixed state in our model, it is not clear 
how general this conclusion is, i.e. whether it is valid be- 
yond the BdG theory. We think that within the Fermi 
liquid scenario both are viable options at this point and 
further work will be needed to distinguish between them. 
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* timchen@mit.edu 

[1] N. Doiron-Leyraud et al, Nature 447. 565 (2007). 

[2] A. Bangura et al, Phys. Rev. Lett. 100, 047004 (2008). 

[3] E.A. Yelland et al, Phys. Rev. Lett. 100, 047003 (2008). 

[4] C. Jaudet et al, Phys. Rev. Lett. 100 187005 (2008). 

[5] S.E. Sebastian et al, Nature 454, 200 (2008). 

[6] S. Chakravarty and H.-Y. Kee, PNAS 105, 8835 (2008). 



[7] A.J. Millis and M. Norman, Phys. Rev. B 76, 220503(R) 
(2007). 

[8] J. Dimov et al, Phys. Rev. B 78, 134529 (2008). 

[9] D. Podolsky and K-Y Kee, arkv 0806.0005. 
[10] D. LeBoeuf et al, Nature 450, 533 (2007). 
[11] J.M. Harris, N.P. Ong, and Y.F. Yan, Phys. Rev. Lett. 

71, 1455 (1993). 
[12] S.H. Pan et al, Phys. Rev. Lett. 85, 1536 (2000). 
[13] Lu Li et al, J. Magn. Magn. Mate 310, 460 (2007). 
[14] T. Janssen et al, Phy. Rev. B 57, 11698 (1998). 
[15] J. Wosnitza et al, Phy. Rev. B 62, 11973 (2000). 
[16] K. Maki, Phys. Rev. B 44, 2861 (1991). 
[17] M.J. Stephen, Phys. Rev. B 45, 5481 (1992). 
[18] T.C. Hsu, B.J. Marston, and I. Affleck, Phys. Rev. B 43, 
2866 (1991). 

[19] S. Chakravarty et al, Phys. Rev. B 63, 094503 (2001). 
[20] PA. Lee and D.S. Fisher, Phys. Rev. Lett. 47, 882 
(1981). 

[21] S.E. Sebastian et al, unpublished. 



